Critical rotation of an annular superfluid Bose gas 
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We analyze the excitation spectrum of a superfluid Bose-Einstein condensate rotating in a ring 
trap. We identify two important branches of the spectrum related to outer and inner edge surface 
modes that lead to the instability of the superfluid. Depending on the initial circulation of the 
annular condensate, either the outer or the inner modes become first unstable. This instability 
is crucially related to the superfluid nature of the rotating gas. In particular we point out the 
existence of a maximal circulation above which the superflow decays spontaneously, which cannot 
be explained by invoking the average speed of sound. 
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After the pioneering work on persistent flow in he- 
lium [H, recent experimental success at producing circu- 
lating superfluid flow of Bose gases in annular traps [3- 
0] has focused interest on the issue of dissipation of this 
macroscopic quantum state. In a superfluid this question 
is crucially related to the existence of a critical velocity Vc 
above which excitations are generated. The critical veloc- 
ity is determined by the Landau criterion [Hj . Dissipation 
occurs for a fluid velocity larger than v^. Symmetrically 
a defect moving above Vc generates excitations in a su- 
perfluid at rest. This has been evidenced experimentally 
in trapped Bose gases @]. 

In a homogeneous weakly interacting Bose gas Vc is 
equal to the speed of sound ■ This is no longer true if 
the system is inhomogeneous. For example, in an infinite 
cylindrically symmetric tube with transverse harmonic 
confinement, the critical velocity is lower than the speed 
of sound Q . In such a geometry, the first modes excited 
at the critical velocity have been shown to be surface 
modes propagating along the tube and localized sym- 
metrically all around the edge of the cylinder. 

In order to study superfluidity experimentally, it is nat- 
ural to bind this system and investigate the stability of 
a persistent flow in a ring geometry. A crucial difference 
between a tube and a ring is the presence of a centrifu- 
gal force arising from the non-Galilean nature of rota- 
tion. Moreover, the curvature of the annulus makes the 
inner and the outer surfaces of the fluid no longer equiv- 
alent ^ (see Fig. [J). 

The ring geometry has recently attracted a lot of in- 
terest. Many annular traps have been proposed [TT| - [l3| 
and realized 0-i, Q [il- Studies of the superfluidity 
include the observation of a persistent current 0, the 
effect of a weak link Il7l |. and the observation of 

a stepwise dissipation of the circulation [3| • The ground 
state of the system in the presence of rotation has been 
determined theoretically [18l - l2C)l | . Phase fluctuations in a 
ring trap have also been investigated [2l|. However, the 
determination of the critical angular velocity in a ring 




FIG. 1. (Color onUne) Sketch of the system: A Bose gas 
flowing with a circulation I (dashed line) and a Thomas-Fermi 
width 27? is held in a ring trap of radius ro. A defect rotating 
counterflow at angular velocity Q, above the critical velocity 
will induce dissipation. 



is still an open question and is highly relevant to recent 



is still an open questic 
experiments 0, 0, llBl- 



* |helene.perrin@univ-parisl3.fr| 



In this Rapid Communication we determine the critical 
angular velocity in the sense of the Landau criterion for a 
Bose gas trapped in a ring. We compute the Bogoliubov 
excitation spectrum both for an initially non-rotating gas 
in the ground state, and for an initially circulating sta- 
tionary state. We show that the critical velocity is gov- 
erned by surface modes, like in the case of an infinite 
tube. However, we find that there are now two distinct 
nondegenerate families of surface excitations propagat- 
ing either at the inner or at the outer edge. The lowest 
of these two branches gives the critical angular veloc- 
ity. Our numerical calculations predict the existence of 
a maximal circulation, above which the system becomes 
unstable, any static perturbation giving rise to dissipa- 
tion of the flow. We give a simple interpretation of all 
these features by extending the surface mode model [i| to 
the ring geometry. Our model is in good agreement with 
the numerical calculations even for an initially circulating 
state. 

We consider a condensate confined in a ring-shaped 
trap (see Fig. [1]) described by the Gross-Pitaevskii equa- 
tion at zero temperature. For simplicity, we reduce the 
problem to two dimensions in the plane of the ring. The 
two-dimensional (2D) ring still allows one to identify the 



2 



inner and outer edges and evidenee their respective roles. 
The trapping annular potential is written as a harmonic 
potential of frequency centered at a radius po- In the 
following, we use the associated scales for energy (hiOr), 
time (a;~^) and length = \Jh/{Mujr), where M is the 
atomic mass. The 2D Gross-Pitaevskii equation reads 



(1) 



where ip = xp{r,0,t) is normalized to unity, A ~ d'^ + 
dr/r + dg/r^ is the Laplacian in polar coordinates (r, 9), 
To = po/or is the dimensionlcss ring radius, and g is 
the dimensionless interaction constant in two dimen- 
sions [28| . 

Using the rotational invariance of Eq. ([T]), we consider 
solutions of the form: 

^(r, 9, t) - e-^(^*-^«) [Mr) + HL{r, 0, t)] , (2) 
where 

<5^i(r,0,t) =^l(r)e-'("*-™'') +«i(r)*e'("**-™'".(3) 

The stationary solution ')pi{r) is a state of circulation £ 
and chemical potential jj,, which depends on i, and Sipf^ 
is a small perturbation, parametrized by £ and m. '(/'^(r) 
is not necessarily the ground state of the system but can 
be realized experimentally using phase imprinting d, 0] . 
We label as i? = \/2jl the half width of the radial density 
distribution in the Thomas-Fermi approximation. 

Linearizing Eq. (JlJ, we solve the Bogoliubov-de 
Gennes equations for and v^{r) [231. We get 

real frequencies with a dispersion relation uj = uJi{m) 
for each initial circulation £. The lowest branch of the 
spectrum allows us to compute the critical angular ve- 
locity ricii) = mium [uJe{m)/ \m\] for a given circulation 
£. 

Figure [2] shows a typical spectrum obtained from 
Eq. ^ for a non-circulating initial state {£ — 0). At small 
TO values, the lowest branch is linear, iooijn) = mQs, and 
can be associated to rotating soundlike waves with an- 
gular velocity fig. At larger to values, this branch ex- 
hibits a small negative curvature that makes the criti- 
cal angular velocity smaller than the angular speed of 
sound [r2c(0) < ^s]- This is not surprising as it is al- 
ready the case in a linear geometry for an inhomogeneous 
gas |@, d^l- As TO increases, the radial profile of the asso- 
ciated density perturbation Sp^{r) of the lowest energy 
mode, where 6pi,{r) = 2Re[Mr)* {uL{r) + vl{r)*)] , 
is more and more localized on the outer radius [Fig.SJJb)]. 
We thus expect that the mode corresponding to the crit- 
ical angular velocity will be correctly described by a sur- 
face mode model. 

Following the approach of Ref. Q, we find the criti- 
cal velocity for a family of modes lying on the edge of a 
condensate. Locally, the surface can be considered as a 
plane, and a surface excitation with wave vector k parallel 
to this plane gives a critical linear velocity Vc — V^p^^^ 
in our dimensionless units, for the critical wave vector 
fee ^ 0.89 X 72/^1/6 d]. In our ring-shaped geometry 
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FIG. 2. (Color online) Excitation spectrum obtained from 
Eq. IIJ (see text for details), with £ = 0,ro = 12 and g = 9000. 
Only the lowest eleven branches are shown. The solid lines 
are a guide to the eye to distinguish between the different 
branches. A symmetric spectrum also exists for negative val- 
ues of m. Solid lines: relation dispersion in the surface mode 
model for the inner (upper blue line) and outer (lower red 
line) modes. 



we identify the critical angular velocity as flc = Vc/fe, 
where is the radius at which the excitation is local- 
ized. Within the surface mode model, the critical angular 
velocity is then: 



y2/ii/6 



(4) 



kcTe, where 



and corresponds to a critical excitation to, 

— tq + R (respectively tq — R) for an excitation lying 
on the outer (respectively inner) edge of the condensate. 
The dispersion relations of the inner and outer edge sur- 
face modes are plotted as solid lines in Fig. [2j For an 
initial state ^ = 0, the inner mode belongs to a higher 
branch and thus does not determine the critical velocity. 

Figure [3] shows a comparison of the full numerical cal- 
culation of the critical velocity with the surface mode 
model. The model of Eq. ^ is in good agreement with 
the numerical calculation and can thus be used to get 
an estimation of the critical velocity. We note that the 
agreement with the numerical calculation is better for 
larger p, as shown in the inset of Fig. [31 since the exci- 
tation is sufficiently localized on the surface @. 

We now turn to the more complex situation of an ini- 
tial state with given circulation £, as obtained experimen- 
tally in recent experiments 0] ■ Figure [H^a) shows the 
lower branch of the excitation spectrum for initial states 
with increasing circulation. As expected, for ^ > 0, the 
spectrum becomes asymmetric, excitations with to < 
propagating against the superflow having lower energies 
than those with m > 0, [Fig. [4l^a), inset]. One striking 
feature of this spectrum is that for a sufficiently large 
£, the critical velocity is associated with a branch that 
crosses the outer edge surface mode branch. The radial 
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FIG. 3. (Color online) Inverse of the critical angular veloc- 
ity as a function of the ring radius ro, at fixed chemical po- 
tential [30[, as obtained from the full numerical calculation 
(dots). Inset: critical angular velocity as a function of fi for 
ro = 40. In both graphs, the solid line is the model of Eq. Q 
with re = ro -I- R. 
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profile (5/9™ (r) shows that these modes are located at the 
inner edge of the ring [see Fig. |3Jb)]. Hence, depending 
on the initial circulation, the most probable mechanism 
for dissipation implies either outer or inner modes. In- 
terestingly, the phase profile of the perturbation displays 
antivortex patterns on the outer edge for outer modes 
and vortex patterns on the inner edge for inner modes. 
This supports the idea that these modes are precursors 
of (anti) vortex nucleation. The total density [i/'l^ with 
a fraction of 10% in the excited mode is displayed in 
Fig. lU^c) for the inner mode, and Fig. IH^d) for the outer 
mode. In both cases, vortices (antivortices) are located 
on the inner (outer) edge of the gas where the density 
vanishes. 

This result can be understood by extending the surface 
mode model to an initial state with circulation. It is 
important to remark that due to the superfluid nature of 
the condensate flow, the local velocity at the inner edge 
is larger than the one at the outer edge. In the frame 
corotating with one of these edges, where the condensate 
surface is at rest, the result of the surface mode model 
still holds. In the laboratory frame, the critical angular 
velocity is then shifted by the angular velocity of the 
corotating frame: ilr = i/r'^, which depends on the edge 
considered, and may be written 



(5) 



The result of Eq. ([5]) is the sum of two contributions: the 
first one arises from the surface mode model, whereas 
the second one is due to the superfluid rotation of the 
condensate itself. These two terms depend with differ- 
ent power laws on the excited mode radius rg and this 
feature explains the transition between inner and outer 
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FIG. 4. (Color online) (a) Lower branch of the excitation 
spectrum for m < 0, ro — 12, g = 9000 and initial states 
with increasing circulation, ^ = 0, 6, 11, 12, 13, 15. Inset: full 
spectrum for £ — 11 (black solid line) and surface mode model 
relation dispersion for outer mode (red dashed line) or inner 
mode (blue dotted line), (b) Radial density profiles of the 
condensate (black solid line), sound-like mode {£ = 0, m = 
— 1, green dash-dotted line, scaled x3 for clarity), inner edge 
surface mode {£ — 15, m = —14, blue dotted line) and outer 
edge surface mode {£ = 0, m = —34, red dashed line), (c-d) 
Corresponding plot of the atomic density with a population 
of 10% in the inner mode (c) and the outer mode (d). The 
size of both images is 40 x 40 in units of a^. At the periphery 
of the gas, 14 vortices appear in (c) and 34 anti-vortices in 
id). 



edge surface excitations observed in Fig.UJ^a) (see inset). 

Figure [5] shows the critical angular velocity as a func- 
tion of the initial state circulation £. The curves exhibit a 
piecewise linear dependence on the circulation. The two 
slopes correspond, respectively, to modes lying on the 
outer or inner edge. The critical angular velocity and 
these slopes are compared to the surface mode model. 
The agreement is good except for small radii where the 
centrifugal term l/r in the Laplacian plays an impor- 
tant role, especially in the case of the inner mode. At 
large values of tq the difference between inner and outer 
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FIG. 5. (Color online) Critical angular velocity versus ini- 
tial state circulation, as obtained from the numerical solu- 
tion (dots), for ro — 20 and g = 15 000. Inset: slope of the 
critical velocity for inner (blue squares) and outer (red dots) 
modes versus the ring radius, at fixed chemical potential [s^l • 
For both graphs, the lines are the predictions of Eq. ((SJ, for 
Te = ro + R (solid line) and re = tq — R (dashed line). 



modes becomes less pronounced as the system resembles 
more and more an infinite tube, where these modes are 
degenerate. 

To interpret the data we define two thresholds. The 
first threshold £i = /^^/^(rg - R^)/{y/2ro) is obtained 
when Eq. ([5]) gives the same value for rg = rg ± R. It 
corresponds to the frontier between the domains where ei- 
ther the inner or the outer modes govern dissipation. We 
find that above a second threshold £2 = ^^^^J'^^^{'ro — R), 
for which the critical angular velocity vanishes, the sys- 
tem is unstable, in the sense that the computed spectra 
contain negative energies. Any static perturbation of the 
system thus triggers dissipation and circulation becomes 
highly unstable. It is important to point out that the 
value of £2 is related to the velocity of the surface mode, 
not to the speed of sound. 

From a practical point of view, this sets an upper limit 
on the circulation that can be imprinted for a given ring 
radius and chemical potential. Therefore, even if the ring 
geometry is adapted to study the persistent flow of a su- 
perfluid, it can bear only a limited amount of circulation. 

Within this model we can make predictions for the 
superflow stability in recent experiments. In the presence 
of static defects, the superflow wiU be unstable for £ > 
£2- We compute this maximum allowed circulation for a 
stable flow with the experimental parameters of Ref. Q , 
namely, ro = 10 and n = 10.9, and find £2 — 11- We may 
wonder if our 2D model applies to three-dimensional (3D) 
experiments. The surface mode model depends only on 
the force at the surface For a harmonically trapped 
BEC, this force is fully determined by the Thomas-Fermi 



radius and the oscillation frequency. A 2D model is thus 
expected to compare well to a 3D experiment with the 
same Thomas- Fermi radius (or chemical potential). 

Our results enlighten the recent work of Ref. [l^ where 
a dynamical simulation of a circulating annular Bose- 
Einstcin condensate in the presence of a static weak link 
shows a dissipation mechanism based on two critical bar- 
rier heights, associated to a vortex-antivortex annihila- 
tion. Our work shows that indeed a static defect can 
induce dissipation by first coupling to an inner edge sur- 
face mode, allowing a vortex to nucleate, as shown in 
Fig. m^c). The excitation of the outer edge surface mode, 
implying anti-vortices located all around the outer edge 
[see Fig. |4l[d)], requires a stronger excitation. 

In conclusion we have computed the Bogoliubov spec- 
trum of circulating annular Bose gases and obtained ana- 
lytical expressions for the critical angular velocity based 
on a surface mode model. We have pointed out the role of 
inner and outer modes in the determination of the critical 
angular velocity. We have discussed the implications of 
these results to explain the dissipation of a persistent flow 
and have shown that the circulation is unstable above a 
given threshold. In fact, the Landau argument Q is quite 
subtle in a ring geometry. Indeed, as shown in this paper, 
we find different results for the motion of a defect in a 
superfluid at rest or for a circulating superfluid flowing 
through a static defect. In the former case outer edge 
surface modes are always excited first, while in the lat- 
ter case the inner edge surface modes dominate. This 
suggests that in an annular geometry the notion of local 
speed of sound, associated with modes centered at the 
peak density, is not the most relevant to discuss super- 
fluidity. 

Further work will include a numerical study of the dy- 
namics of an annular Bose gas in the presence of static 
and rotating defects to further investigate the critical 
velocity. In particular, we expect that engineering the 
shape of a perturbation may help to selectively excite 
only the critical mode and thus more precisely control 
the system. An interesting point would be to examine 
the possibility of inducing a circulation by selectively ex- 
citing this mode. Symmetrically, while it is clear from 
our results that dissipation of a flow with an initial cir- 
culation £ > £2 flrst occurs through the appearance of 
vortices at the inner edge, the subsequent evolution is an 
interesting open question which will require the simula- 
tion of the full dynamics. 
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